Mon. Not. R. Astron. Soc. 000, [TpS] (2011) Printed 24 February 2013 (MN I^TeX style file v2.2) 



The JCMT Nearby Galaxies Legacy Survey. VII. Ha 
imaging and massive star formation properties 
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ABSTRACT 

We present Ha fluxes, star formation rates (SFRs) and equivalent widths (EWs) for 
a sample of 156 nearby galaxies observed in the ^^CO J = 3 — 2 line as part of the 
James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey. These are derived 
from images and values in the literature and from new Ha images for 72 galaxies which 
we publish here. We describe the sample, observations and procedures to extract the 
Hof fluxes and related quantities. We discuss the SFR properties of our sample and 
confirm the well-known correlation with galaxy luminosity, albeit with high dispersion. 
Our SFRs range from 0.1 to llMQyr"^ with a median SFR value for the complete 
sample of 0.2 M© yr~^. This median values is somewhat lower than similar published 
measurements, which we attribute, in part, to our sample being Hl-selected and, thus, 
not biased towards high SFRs as has frequently been the case in previous studies. 
Additionally, we calculate internal absorptions for the Ha line, A{Yia)^ which are lower 
than many of those used in previous studies. Our derived EWs, which range from 1 
to 880 A with a median value of 27 A, show little dependence with luminosity but rise 
by a factor of five from early- to late-type galaxies. This paper is the first in a series 
aimed at comparing SFRs obtained from Ha imaging of galaxies with information 
derived from other tracers of star formation and atomic and molecular gas. 

Key words: galaxies: star formation - galaxies: ISM - infrared: galaxies - submil- 
limetre: galaxies. 



1 INTRODUCTION 

Along with galaxy interactions and the presence of active 
nuclei, star formation efficiency is one of the most impor- 
tant factors that controls a galaxy's structure and evolution. 
Although the subject has been extensively researched both 



in the local Universe (e.g.,|Kennicutt fc Kent 



cutt"2008') and in the more distant Universe (e.g., [Madau 



et al |[T995l |Alonso et al.||1999| [James et al 



1983| Gallego 



[20041 IKennT 



et al. 1996; Hopkins, Connolly fc Szalay|200^|Lanzetta et al. 



2002 



Nagamine et al. 2006), many fundamental questions 



in galactic and extragalactic research remain. 

In the local Universe, many surveys have been carried 
out in different wavelengths (see, for example, jKennicutt 



2008 for an extensive compilation of literature data) to probe 



the star formation rate. Among these, the Ha line is one of 
the most widely used star formation rate (SFR) tracers, both 
photometrically and spectroscopically ( [Kennicutt| |1998| . In 
contrast to, e.g., the Far Ultra- Violet (FUV), Ha lies in the 
optical region for nearby galaxies, thus allowing easy ground- 
based observations. In addition, the spatial resolution which 
can be obtained with Ha is superior to what is normally 
obtained at other wavelengths such as the Infrared (IR), 
which are also more complicated to observe with ground- 
based telescopes. Studies carried out throughout almost two 
decades have proved that the Ha line is a solid tracer of star 
formation, especially if the measurements are complemented 
by the use of Far Infrared (FIR) data (e.g., ^Calzetti et al.[ 
2007[|2010|[Leroy et al.||2012[ ). 
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Recent studies (see Wilson et al. 2009 and references 
therein) state that the warm, dense molecular hydrogen is 
better correlated with star formation than the cold, more 
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dispersed gas. Due to the lack of a permanent dipole mo- 
ment, all the ground-based observations aimed to detect 
molecular hydrogen are made using alternative tracers. Al- 



though not the only method at hand (Gautier et al.||1976 



Neufeld & Yuan 2008), the spectral lines produced by the 



different rotational transitions of the CO molecule stand as 
the best known and most widely used tracers of molecular 
hydrogen. However, the vast majority of surveys aiming to 
study the molecular gas in nearby galaxies have focused on 
the lower transitions [J — 1 — 0, J — 2 — 1) of CO, which 
do not discriminate between cold and warm gas. 

The James Clerk Maxweh Telescope ( JCMT) on Mauna 
Kea, Hawai'i is currently hosting a series of Legacy Surveys, 
including the Nearby Galaxies Legacy Survey (NGLS) which 
aims to map 156 galaxies in the CO J — 3 — 2 spectral line 
and in the 450 and 850 /xm continuum (C. Wilson et al., 
in preparation). All the CO J — 3 — 2 observations have 
been completed and the continuum observations have begun. 
Further information about the scope and science goals of 



this survey can be found in Wilson et al. ( 2009 ) as well as in 



the subsequent series of papers which make use of the CO 
J = 3 - 2 d ata ([Warren et al.pOlOllBendo et al.|2010j [Irwin 



et al.|2011[ [Wilson et al.|2011[[Sanchez-Gallego et al.J2011 ) 
Because of the physical conditions needed to activate 
this transition, the J = 3 — 2 line is perfect to trace and 
study the warm molecular gas which is most likely to form 
stars. This directly inspires the interest in observing this 
same sample using a well proven tracer of star formation 
such as Ha. 

This paper is organized as follows. Section [2] introduces 
the sample and the criteria used for its selection; Sections 
|3] and |4] describe the observations and the origin of the lit- 
erature data used in this work. Section |5] explains in detail 
the complete process of data analysis, from the most basic 
reduction to the continuum subtraction, flux measurements 
and SFR and EW calculations. In Section |6] we analyze the 
SFR and equivalent width (EW) results obtained. They are 
compared with previous data from the literature in Section 
|7| Finally, Section [S] presents the main conclusions and out- 
lines the future work we plan to do with these data. 



2 SAMPLE SELECTION 

This work presents Ha photometry for a total of 156 nearby 
galaxies. The sample of galaxies has been mostly Hi-flux 
selected in order to avoid a SFR-driven selection (as would 
have happened using a FIR flux criterion) while ensuring 
that the galaxies have a rich interstellar medium. All our 
galaxies are contained within a range between 2 and 25 Mpc 
in distance {v < 1875 km s~^) and were selected to cover 
the complete morphological spectrum, from early- type spi- 
rals and ellipticals to late- type spirals and irregular galaxies. 
Both field and Virgo cluster galaxies are well represented 
in each morphological bin. Most of the Virgo galaxies have 
D25 < 4arcmin, which allowed us to use the jiggle map 
method during the CO J — 3 — 2 observations (Warren et al. 



2010) carried out with the JCMT. C. Wilson et al. (2011, in 



prep.) give a more detailed explanation of the criteria and 
parameters used to construct the sample. 

The complete sample has been compiled from three dif- 
ferent subsets: (1) 48 galaxies taken from the Spitzer In- 



frared Nearby Galaxies Survey (SINGS; Kennicutt et al. 
2003 ); (2) 36 galaxies chosen from a Hi flux-limited sam- 
ple of galaxies from the Virgo cluster; and (3) 72 H i flux 
limited field galaxies. The fact that a significant part of our 
list of targets matches the SINGS sample ensures that a 
large amount of multi- wavelength data is available for those 
galaxies, from the infrared imaging obtained by the SINGS 



team to Hi data (for example, from THINGS; Walter et al 



2008) and UV imaging from GALEX. However, SINGS used 



FIR emission to constrain its sample, which leads to a cer- 
tain bias towards spiral, star- forming galaxies. This bias also 
appears in our sample, as can be seen in Figure [l] which 
shows the distribution of the galaxies of our sample across 
the Hubble sequence and the contribution of each subsam- 
ple. This bias will also become apparent when we discuss 
the SFRs (Section p|. The morphological breakdown is as 
follows: 4 elliptical galaxies (T G [—6,-4)), 20 lenticular 
galaxies (T G [-4, 0)), 50 early spiral galaxies (T G [0, 4.5)), 
49 late spiral galaxies (T G [4.5, 9.5)) and 33 irregular galax- 
ies (T G [9.5,10]). We make a distinction between small 
(D25 < 4arcmin) and large galaxies in the SINGS subsam- 
ple. The sample is given in Table 1. 



3 NEW Ha OBSERVATIONS 

We have imaged 72 out of the 156 sample galaxies in the 
Ha line. While a small subset of these galaxies had Ha data 
available in the literature, we decided that the images were 
either not deep or not detailed enough to serve our purpose 
and we reobserved them. 

For these observations, we made extensive use of sev- 
eral of the telescopes located on the Roque de los Muchachos 
Observatory (La Palma, Canary Islands). In particular, we 
used: (1) the Wide Field Camera (WFC) mounted on the 
2.5 meter Isaac Newton Telescope (INT), (2) the Auxiliary 
port Camera (ACAM) and the auxport camera (AUX) on 
the 4.2 meter Wilham Herschel Telescope (WHT), and (3) 
the Andalucia Faint Object Spectrograph and Camera (AL- 
FOSC) on the 2.5 meter Nordic Optical Telescope (NOT). 
Additional images were taken using the MDM8K mosaic 
camera on the 2.4 meter Hiltner telescope at the Michigan- 
Dartmouth-MIT Observatory (Kitt Peak). Table 1 summa- 
rizes the dates, observing conditions, exposure times and 
filters used for each galaxy. 

Although working with different telescopes and instru- 
ments, we kept a systematic observation strategy. Each 
galaxy was observed through a narrow-band filter redshifted 
to the position of the Ha line for each galaxy. The narrow 
filters used had a width in the range 15 — 50 A with trans- 
missions always above 50% and, generally, higher than 70%. 
We ensured that the Ha line always fell in the region of max- 
imum transmission of the filter. For each target we obtained 
at least three Ha images in order to remove cosmic ray hits 
through the combination of the different images. Standard 
exposure times were 3 x 300 s. 

The continuum subtraction was done using an R broad- 
band image. Due to the low recession velocity of our galaxies, 
no other filters were necessary to obtain the continuum con- 
tribution in the Ha region. We typically obtained 3 x 180 s 
continuum exposures which were combined to produce a sin- 
gle i?-band image. 
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Sky flat frames were taken during both the evening and 
morning twihghts. In some cases this was not possible due to 
adverse atmospheric conditions during one or both twilights. 
In some of these cases, dome flats were taken. A series of 
bias exposures was obtained at the beginning of the night 
and again at the end. Typically the bias was constant and 
uniform throughout the night. 

The photometric calibration was carried out by observ- 
ing at least one spectrophotometric standard star at the 



beginning of the night and at the end. Landolt ( 1992 ) pho- 



tometric stars where used for the R-hand calibration. While 
most of the images were obtained during photometric nights, 
some galaxies were observed in clear but not perfect condi- 
tions. In those cases we did not calibrate the images. Instead, 
we reobserved the galaxies during photometric nights, ob- 
taining short exposures (typically 1/4 the original exposure 
times) which we used to calibrate the previous images. 



4 LITERATURE DATA 

The remaining 84 galaxies in our list had previously been 
observed, and photometrically calibrated Ha imaging was 
available in the literature. These include: 25 galaxies from 
the Ha Galaxy Survey (HaGS; James et al.|2004); 33 galax- 



ies from SINGS (Kennicutt et al. |2003); 4 galaxies from 



Knapen et al. ( 2004 ); 4 galaxies from|Kennicutt et al.,(|2008 ); 



and 18 galaxies from the GoldMine project ([Gavazzi et al, 
2002[ |2006| [Boselh et al.||2002t 



For these galaxies, the headers of the FITS files have 
been modified to match the global standards of this work 



(see Section 5.4) and the images have been divided by the 
exposure time to get an effective exposure time (texp) equal 
to 1 second. If necessary, the astrometric solution in the 
images was improved. 

When possible (i.e., when both the original Ha and 
the continuum images were available) we did not use the 
Ha continuum-subtracted image from the literature, but we 
performed our own continuum subtraction as explained in 
Section [5^ We used the final images to measure our own 
flux values and derived properties, as described in Sections 
15.515.81 



5 DATA REDUCTION 

5.1 Basic reduction 

The data reduction is quite straightforward and we refer 
to the corresponding procedures described in |Knapen et al.| 
(2004). The reduction was performed using mainly packages 
and tasks from IRAFJj Again, slightly different procedures 
were carried out depending on the instrument used to ac- 
quire the images. The main steps were, however, common 
to all the images regardless of their origin. The basic steps 
were as follows. 



a) A median bias frame was produced combining all valid 
bias exposures. This median bias was then subtracted from 
each of the flat field and science images. Some images from 
the literature were bias-corrected using a constant value 

^ http://iraf.noao.edu/ 

© 2011 RAS, MNRAS 000, [TpS] 



(namely the mean value of the median bias frame) rather 
than the master bias image (see Knapen et al.__2004j . This 
is no problem as the structure of the individual bias frames 
was negligible in those cases. 

b) The flat field correction was applied using an adequate 
median sky flat obtained from the flat frames for each fil- 
ter. Evening and morning sky flats were analyzed separately 
although, in the majority of cases, no important differences 
were found and all the flat frames were combined. 

c) The images were aligned, their background was sub- 
tracted and they were weighted by their exposure time and 
then combined by filter using the median. A sigclip rejection 
algorithm was applied to remove pixels with values above or 
below 3a. This ensures that most of the cosmic rays were 
successfully removed. 

d) Although in most cases the images had a good astro- 
metric calibration, a new astrometric solution was applied 
to each image using the Guide Star Catalogue II (GSC-II; 



Lasker et al. 2008). The final precision is better than one 



5.2 Continuum subtraction 

Once both the Ha on-line and the continuum images are 
completely reduced we must subtract the contribution of the 
continuum to the Ha flux. The great majority of our con- 
tinuum images were obtained using i?-band filters and the 
rest of this subsection will deal with that situation. How- 
ever, for a small number of cases (all from the literature) 
the continuum image is a narrow filter of width similar to 
that of the Ha filter but shifted in wavelength so that only 
the continuum and not the line is measured. In these cases 
the procedure to subtract the contribution of the continuum 
is similar to when using i?-band, but much more straightfor- 
ward as the scale factor between the Ha and the continuum 
images is always ^ 1. 

The method used to calculate the scale factor of the 



continuum image is well explained in Knapen et al. (2004) 



and we refer to that paper for further details. First, if there 
is an important difference between the full widths at half 
maximum (FWHMs) of the Ha and continuum images, we 
smooth the images to the worst FWHM using a Gaussian 
convolution. Then, we plot the intensity of each pixel in 
the Ha image, in counts per second, versus the intensity 
of the same pixel in the i?-band image. In general, a cer- 
tain percentage of the brightest pixels (usually around 3%) 
is removed to avoid the contribution of field stars. The re- 
sulting plot shows a strong linear correlation which can be 
fitted using a least squares regression line, the slope of which 
gives the scale factor between the two images. Although this 
method has been found to produce consistently good results, 
some galaxies (especially those with very low emission) re- 
quire further processing to obtain the optimum scale factor 
value. Thus, a final small correction by hand was often made 
to ensure that the continuum subtraction was optimal by 
looking at different regions of the galaxy and ensuring that 
no area was oversubtr acted. While it is not always trivial 
to find the most accurate scale factor, small differences in 
the value of this parameter do not affect the derived values 
(flux, SFR, EW) in a significant way. 

It must be taken into account that the continuum im- 
age used is contaminated by the flux of the Ha line. The 
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contribution of the line has previously been estimated to be 



3% of the total integrated flux in the i?-band (James et al 
2004). This effect does not occur when the filter used for the 



continuum subtraction is a narrow one, as in that case the 
Ka line is not present at all. Integrating the transmission 
curves of the filters used, we calculated how much the Ha 
line contributed to the continuum flux. We found that, for 
most sets of filters, the correction was in the range 2-4% 
with some extreme cases as high as 6%. The calibration was 
modified according to the value found for each galaxy. Table 
2 shows the correction applied. 

5.3 Photometric calibration 

The next stage in the reduction process was to calibrate the 
science frames photometrically. The i?-band calibration was 



done using Landolt (1992) standard stars and taking the 
airmass of the objects into account. For the calibration of 



the Ha images, following Knapen (1992), we considered the 
flux density, /a, given by 



log /a = -0.4mag;, 

-2 X-1 



.42, 



with /a in units of ergs~^ cm~^ A~^. To calculate the zero 
point factor we performed aperture photometry of the spec- 
trophotometric standard stars selected. For these stars, the 
value of mag^ is tabulated every few angstroms. All the stars 
used are from either (Stone (1977), Massey et al. (1988) or 
Oke ( 1990 ). If X is the total flux of the standard star in units 
of counts s~^, the zero point factor can be written as 

AA/a 



ZP: 



X 



where AA corresponds to the width of the Ha filter used. ZP 
has units of ergcm"^ count"^. The flux for a certain region 
can be calculated as 

logioF = ZP + logio (counts), 

in units of ergcm~^s~^. In this expression we have taken 
into account that the effective exposure time of our images 
is one second. 



5.4 Final image preparation 

The final images for each galaxy (Ha without continuum 
subtraction, i?-band and continuum-subtracted Ha) were 
rotated North-up and trimmed. The FITS headers of all 
images (including those obtained from the literature) have 
been modified in the following way: first, a custom header 
has been inserted including the relevant data (total expo- 
sure time, zero point, flux and luminosity per count, etc); 
second, the original header has been preserved after our cus- 
tom header. The units in the final images are always counts 
per pixel (i.e., the exposure time has been set to one sec- 
ond). The original exposure time for each frame is stored in 
the header of the image under the field darktime. 

When working with literature data, the astrometric so- 
lution provided with most of the images was well below the 
FWHM of the frames and no correction was made. In some 
cases, however, the astrometric solution was wrong or nonex- 
istent. Those images were astrometrically calibrated using 
the Guide Star Catalogue (Lasker et al. 1990) which pro- 
vides an accuracy better than one arcsec. 



5.5 Flux measurement 

In order to determine the Ha flux without the contribution 
of field stars, we masked all the stars near the galaxy while 
taking special care not to remove any H II regions. We used 
the masking utility integrated in Gaia ( Draper et al.|[2009 ) 
which substitutes the region to mask based on an estimate 
of the values of the surrounding pixels. We applied a third 
order surface fit, which in all cases produced an excellent 
result. The stars were masked in the continuum band and, 
when necessary, also in the Ha continuum-subtracted image. 
To calculate the integrated Ha flux of the galaxies we 
used the task ellipse in IRAF. On the continuum sub- 
tracted frame we ran ELLIPSE to calculate the integrated flux 
for a series of concentric ellipses centred on the nucleus of 
the galaxy and with ellipticities (s) and position angles (PA) 
as obtained from the Third Reference Catalogue of Bright 
Galaxies (RC3; de Vaucouleurs et al.|1 99r). For a number of 
galaxies, especially in the case of irregular ones, these data 
were either not available or not reliable. In these cases circu- 
lar apertures were used. Table 2 includes the values of e and 
PA we used in our measurements. The separation between 
the different ellipses is a function of the pixel scale of the 
image used. Typically we calculated ellipses every 5 arcsec 
although this value was manually modified for special cases 
such as extremely large or small galaxies. The position of 
the nucleus was determined in each case from the contin- 
uum image, except for some irregular galaxies where the 
estimated centroid of the galaxy was used. The determina- 
tion of the optimum semi-major axis of the ellipse was done 



using the same criterion as James et al. (2004): we consider 



that the photometry has converged when the measured flux 
between two consecutive ellipses varies by less than 0.5%. 
This method works efficiently for most of the galaxies in our 
sample. However, for a small number of galaxies (mainly ir- 
regular ones) the resulting aperture was clearly wrong. In 
those cases we used the D25 value from the RC3. The aper- 
ture determined in this way was then used on every available 
band. The fluxes obtained were converted from counts to 
erg cm~^ s~^ (in the case of the Ha band) or to magnitudes 
(case of i?-band) using the previously calculated conversion 
factors. 

Figures [2] to [5] show the final Ha continuum-subtracted 
images. In each case, the elliptical aperture where we mea- 
sured the Ha flux is plotted on the image. 



5.6 [Nil] correction 

Although for the sake of clarity we have only discussed 
Ha fluxes, a number of our Ha images (mainly those ob- 
served using filters wider than 35-40 A) are contaminated 
with emission from the [Nil] lines at AA6548, 6584 A. We 
estimated the contribution of the [Nil] lines using the ex- 
pression from Kenn icutt] ( |2008 ) 

r -0.173 ±0.007 Ms- 
log ([Nii]/Ha) = <^ -(3.903 ±0.137) if Mb > -21 
[ 0.54 if Mb ^ -21 

(1) 
Helmboldt et al. ( 2004 ) provide a similar expression but 
We decided against using this 



using Mr instead of Mb. 

latter reference because all but two of our galaxies have 

measured values of Mb, while many do not have reported 
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Mr. For those cases where both Mb and Mr measurements 
were available we compared the [Nil] /Ha obtained using 
both methods, finding a good agreement. [James et al. ( 2005 ) 



studied the [Nil] /Ha ratio for the HaGS galaxies obtaining 
values which are compatible with ours. 

Even when it is possible to know the relation between 
[N ll] and Ha using Equation [l] we still have to deal with 
the problem of the different transmissions for the Ha and 
the [Nil] lines. Except for the broadest filters, the Ha line 
is located in the region of maximum transmission of the 
filter (see Section [s]) but the [Nil] lines at AA6548,6584A 
lie at a position where the transmission is much reduced 
or even marginal. Solving this problem is not easy as filter 
shapes are often different and difficult to model. Most of the 
filters we used are neither Gaussian nor rectangular but a 
convolution of those two shapes: i.e., the transmission re- 
mains approximately constant for most of the filter width 
and then decays exponentially. To estimate the real trans- 
mission of the [Nil] lines, we took a series of filters with 
different FWHM (from 36 to 90A) for which a transmission 
curve had been accurately measured. We calculated the re- 
lation between the transmission of the filter for the position 
of the Ha and the [Nil] lines, T[Nii]/Ha- Figure [g] shows the 
values obtained. We fitted the points to a sigmoid model ob- 
taining that T[Nii]/Ha = 1 - 2.670 e-^^"^/^^-^^^ Using this 
expression we calculated the observed ratio, ([Nii]/Ha)obs, 
as ^obs = ([N ii]/[Ha]) x T[Nii]/Ha, for each galaxy (see Table 
2). 



5.7 Star formation rates and equivalent widths 

To determine the star formation rate (SFR) from the mea- 
sured fiuxes we applied the following conversion fromJKen^ 



nicutt et al. 



(2009): 



SFR(M0 yr" 



5.5 X 10"^^L(Ha), 



where L(Ha) is the luminosity, calculated as 

L(Ha) [erg s"^] =47ri:>^ (3.086 x lO^'^fF^^, 

with D the distance to the galaxy in Mpc and F^^ the 
absorption-corrected fiux. This determination of the SFR 
conversion factor is based on a "Kroupa" initial mass func- 



tion (Kroupa & Weidner 2003). The SFRs based on this 



estimation are lower than others calculated using older re- 



Kennicutt 



nicutt et al. 



lations (for example, Kennicutt, Tamblyn & Congdon 



[1998 ) by approximately a factor of 1.5 (Ken- 



1994 



2009). The distances used here (see Table 2) 



are the same as those indicated in C. Wilson et al. (2011, in 
prep.). 

To correct for Galactic absorption we used, for each 
galaxy, the total absorption value for the R band, A{R) 
given by NEErJ (see Table 2), which is based on the dust 
maps published by Schlegel, Finkbeiner & Davis (1998). To 



correct for the internal absorption, several solutions have 
been proposed in the literature. A number of studies over 
the years have adopted a constant value A{Ha) = 1.1 mag 
(Kennicutt & Kent 1983), regardless of the morphological 



2 The NASA/IPAC Extragalactic Database; 

|http://ned www. ipac.caltech.edu/ 



type or inclination of the galaxy (see, among other s, [Ken- 
nicutt, Bothun fc Schommer 1984] [Kennicutt 19981 [Jgtmes 



et al.|2004[ ). [Niklas, Klein &; Wielebinski[ ( [1997[ ), however, re- 



ported that the value for A(Ha) might be as low as 0.8 mag 
and James et al. ( 2005| found values as low as 0.4 mag for 
some early-type galaxies. On the other hand, [Helmboldt[ 
et al. (2004) proposed an expression to calculate A (Ha) in 
function of the absolute R magnitude of the galaxy, Mr. 
This expression takes the form 

log A(Ha) = (-0.12 ± 0.048) Mr + (-2.5 ± 0.96). (2) 

We used this equation to calculate the internal absorp- 
tion of the galaxies in our sample. Unfortunately, not all of 
our galaxies have been accurately measured in the i?-band. 
For those galaxies we overcame this problem by selecting 
different values of A{Ha) according with the type of galaxy 
studied. Specifically, we used 



A(Ha) 



r 0.3 
I 1-1 



ff Mb ^ -16, 
otherwise. 



which are in good agreement with the values yielded by 
Equation |2] Our method provides absorption values which 
are somewhat lower than those often used in previous stud- 
ies of the same type. This explains that the mean SFR values 
we obtain here are lower than those frequently found in the 



literature (see Section 6.1). 



The equivalent width (EW) of the Ha line is a direct 
and widely used method to estimate the intensity of the line 
as compared to the emission of the underlying continuum. 
As is usual in the literature, we calculate the EW by form- 
ing a rectangle with the height of the continuum under the 
line and finding the width such that the area of the rect- 
angle matches the area of the spectral line. In practice, we 
calculated the EW of our galaxies applying two different 
methods. First, we used the scale factor obtained during the 
continuum subtraction process, together with the FWHM of 
the Ha narrow-band filter. The equivalent width can then 
be determined as 



EW 



CPSho 



k (CPSc 



CPShc.)/FWHMhc.' 



(3) 



where k is the scale factor and CPSx are the counts per sec- 
ond of the Ha and continuum bands, with CPSna already 
corrected for [Nil] contamination and internal absorption. 
We found that this method yielded very accurate results, 
and our EW measurements are determined mainly using this 
system. Unfortunately, some of the literature data sets in- 
cluded neither the continuum scale factor, nor the original 
Ha image, which is necessary to perform our own continuum 
subtraction. In these cases, we used an alternative method. 
The EW can be simply calculated as the ratio between the 
Ha continuum subtracted fiux and the fiux density in the 
continuum band. As our continuum images were obtaining 
using standard i?-band filters, we converted from R mag- 
nitudes to fiux densities using the parameters which define 
the system (we mostly used Cousins R filters with an effec- 
tive wavelength of 6400 A and a zero point of 3080 Jy; Cox 



2000 ) . As this method relies on physical measurements (with 



the uncertainties associated with every calibration) and not 
all the i?-band filters are identical, we found that the EWs 
determined in this way are less accurate. 

Finally, we used the EW values from the literature 
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where available and perform the calculations described 
above in the rest of the cases. Section [7| compares our EW 
values with those in the literature to study the quality of 
our measurements. Table 2 indicates the final SFR and EW 
measurements for each galaxy and, for the latter, the method 
used for the calculation or the source of the value. 



5.9 Image availability 

With the publication of this article, the Ha and i?-band 
images for all the galaxies in our sample are made publicly 
available through the Strasbourg astronomical Data Center 

(CDS). 



5.8 Uncertainties 

Several factors contribute to the final uncertainties of the 
SFR and EW measurements we present here. First, we esti- 
mated the uncertainty associated with the flat-fielding cor- 
rection. Although this value changes depending on the in- 
strument and even with the specific conditions of the night, 
we find the typical uncertainty to be well in the range 1 — 2% 
with some values as high as 5%. 

Uncertainties in the zero point calibrations during pho- 
tometric nights have typical values of 2%. We used several 
galaxies observed both in photometric and non-photometric 
conditions to estimate uncertainties associated with the lat- 
ter. We find that, even in the worst cases, the uncertainty 
in the calibration is never higher than 5%. Similarly, the 
uncertainty in the flux measurement is typically below 1%, 
with some extreme cases as high as 5% for a few very weak 
galaxies. 

For the uncertainties associated with the distances, we 
found that the values given in NED were frequently un- 
reliable and underestimated the real uncertainty. Except 
for a few galaxies with quality distance measurements, we 
adopted a typical uncertainty of 20% in the estimated dis- 
tance to the galaxy. 

Along with the distance, the main contribution to the 
uncertainty for the vast majority of our galaxies is the con- 
tinuum subtraction. As a galaxy is an extended object there 
is often not a single scale factor applicable to the whole im- 
age. This implies that any method aimed to measure the 
scale factor used in the continuum subtraction is affected 
by a certain amount of subjectivity. To measure the un- 
certainty associated with the removal of the continuum we 
selected 10 galaxies of different morphological types, SFRs 
and EWs. For each galaxy we produced new continuum sub- 
tracted frames using both the first scale factor which led to 
a clearly under-subtracted image and the first producing an 
over-subtracted one. Using these extreme values as a ref- 
erence, we estimated that the uncertainty due to the scale 
factor is in the range from 5 to 15% for most galaxies. Unlike 
some previous studies ( James et al.||2004 | pCennicutt 2008) 
which reported a correlation between EW and the uncer- 
tainty of the scale factor, we find no trend whatsoever in 
this sense. 

Without considering the problem of distance (i.e., tak- 
ing into account only those uncertainties intrinsic to the 
quality of the images) we find that the median uncertainty in 
flux for all our galaxies is ^ 18%. This values is in very good 
agreement with recent surveys using similar techniques. For 
instance, Kennicutt] f 2008 ) obtained a median uncertainty 
(also using the flux measurements) of ^ 16% while [James 
eraLl([2004| reported - 18%. 



6 RESULTS 

Table 2 shows the main results of the photometry of our 
156 galaxies as well as their associated uncertainties. The 
column description is as follows. 

Column (1): the name of the galaxy in the designation 
used in this work. 

Column (2): the morphological class of the galaxy as 
catalogued in NED. 

Column (3): the recession velocity used to calculate the 
Doppler shift of the Ha line, obtained from NED. 

Column (4): the distance to the galaxy considering 
Ho = 70.5kms"^Mpc"\ ^^m = 0.27 and Qvac = 0.73 
and using a model that takes into account the Virgo in- 



fall velocity (Mould et al. 2000). For the galaxies belong- 



ing to the Virgo subsample we have adopted a distance 
of 16.7 ± 1.1 Mpc f Mei et al.||2007t. Additio nally, we used 
the distances from iFreedman et al.l (|2001| for NGC 925, 



NGC2403, NGC 3031, NGC 3198, NGC 3351, NGC 3627 



and NGC 4725; Karachentsev et al. (2002) for NGC 2976, 



NGC 4236, UGC 04305 and UGC 08201; [Tonry et al.| ( |2001 ) 
for NGC 584, NGC 855, NGC 4594, NGC 4736, NGC 4826 



and NGC 5194; .Karachentsev et al. (2004) for NGC 628; 



Dalcanton et al.l (2009) for IC 2574; Macri et al. (2001) for 



NGC 2841; Leonard et al. (p002) for NGC 3184; and "Seth, 



Dalcanton fc de Jong (20051 for NGC 4631. For NGC 3034, 



PGC 023521, UGC 05139, UGC 05336 an d UGC 05423 w e 
used the same distance as NGC 3031 ( [Freedman et al.|2001 ). 

Columns (5): the major and minor axis of the galaxy in 
arcmin. 

Column (6): the B-hand absolute magnitude derived 
from the apparent magnitude in the RC3 catalogue and the 
redshift of each galaxy. The luminosity distance, necessary 
to derive the absolute magnitude, was calculated using a 
model which takes Ho = 70.5, Qm = 0.27 and Qvac = 0.73. 
In a few cases (marked in the table) no B-hand magnitude 
was available in RC3 and we use the closest match from the 
photometric data points gathered in NED for each galaxy. 
For PGC 140287 we used the B-hand magnitude from LEDA 
( |Paturelet"aL][2003l ). 

Column (7): the scale factor applied to the continuum 
image during the continuum subtraction (see Section 5.2). 

Column (8): the estimated contribution of the Ka line 
to the flux measured in the continuum filter. 

Column (9): the total flux of the galaxy in the Ha line 
in units of ergcm"^ s~^. The area used to perform the pho- 
tometry has been selected as indicated in Section [53] This 
flux has not been corrected for any kind of absorption or 
[Nil] contribution. 

Column (10): the total absorption for the galaxy. At, 
in mag. This value comprises the Galactic foreground ab- 
sorption in the R band, A{R), for the position of the galaxy 



taken from Schlegel, Finkbeiner & Davis (1998), and the in- 
ternal absorption, A(Ha) calculated as explained in Section 

EZl 
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Column (11): the estimated ratio ^obs = ([Nii]/Ha)obs 
calculated as explained in Section [5^ 

Column (12): the Ha luminosity of the galaxy, in 
ergs~^. This value has been corrected for absorption using 
At and for the contribution of the [Nil] lines. 

Column (13): the star formation rate for the galaxy in 
units of Mq yr~^ (corrected for absorption and [Nil]). 

Column (14): the equivalent width, in Angstrom, calcu- 
lated as described in Section [5?fl The method used for the 
calculation is indicated in parenthesis meaning (1) that the 
equivalent width is derived from the scale factor used for 
the continuum subtraction as in Equation (l3|; (2) that the 
equivalent width is derived from the fluxes of the images; 
and (3) that we used the value from the literature. In the 
latter case, the reference is always the same as the one listed 
in Table 1 for the corresponding galaxy. 



6.1 Star formation rate measurements 

Figure [7^ shows the correlation between the SFR (in units of 
solar masses per year) and the 5-band absolute magnitude, 
M{B). Our four subsamples (field, small Virgo, and large 
and small SINGS galaxies) are plotted separately in Figure 
|8]to consider the individual properties and possible biases 
of these subsets. 

We find an expected trend, with brighter galaxies host- 
ing a higher amount of star formation. This behavior has 
been previously reported by several authors such as |Tresse| 
fc Maddox| ( |1998 ), who used spectroscopic Ha data from a 
sample of galaxies at z ^ 0.2, and more recently | Knap en] 



&: James] (12009) (their Figure 2), using the sample of the 
HaGS survey (see also Section fTI). Figure 6 in Lee et al. 



(|2009b|), who use d ata of - 300 galaxies from llHUGS^Ken_ 
[nicutt et al.||20Q8| , shows a similar trend for galaxies with 
M{B) ^ -15 (for very faint galaxies, M{B) > -15, the SFR 
seems to drop suddenly to values below 10~^ Mq yr~^). This 



same result is found by Karachentsev & Kaisin (2010) using 
435 galaxies in a volume-limited region of 10 Mpc. The slope 
of our least squares fit is -0.29, somewhat higher than the 
value of -0.41 found by |Lee et al^ (|2Q09b). 

The top right corner in Figure |7| is exclusively occupied 
by galaxies from the SINGS large subsample. The galaxy 
with the highest SFR in our sample (11.31 M© yr~^) is 
NGC 4254 (Messier 99). In the other extreme of the plot we 
find two small irregular galaxies (PGC 23521 and IC3105), 
both with a SFR below 10~^ M© yr~\ 

It is interesting to consider the behavior of the four dif- 
ferent subsamples which compose our sample. Field galaxies 
in the NGLS sample are distributed uniformly along a wide 
range of M(B), which is expected as they have been selected 
based on their H i emission and not on their optical luminos- 
ity. The same can be said for the galaxies contained in the 
SINGS subsamples. The Virgo cluster subsample, however, 
contains galaxies with M(B) well constrained in the range 
[—17, —22] magnitudes. 

We calculate the mediaqj SFR for all the galaxies in 
our sample as 0.211q;^2 ^0 yr~^- This value is consistent 



with those of the different subsamples (0.221qj^2 ^0 y^~'^ 
for the field galaxies, O.IIIqq^ for the SINGS small galaxies, 
and 0.091qq5 M© yr~"^ in the case of Virgo galaxies). Not 
unexpectedly, the SINGS large subsample shows a higher 
median SFR {1.04:'tl)\ll Mq yr~^) as it is mainly populated 
by spiral galaxies (see Figure IT]) with high FIR luminosi- 
ties and, thus, large SFR (see, for example, Calzetti & Ken- 
nicutt 2009). On the other hand, the large uncertainty in 



the SINGS median value shows that this subsample con- 
tains galaxies with a large range of SFRs. The median value 
for the whole sample is lower than those of other studies 
in the literature (see Section l7| . We attribute this differ- 
ence to the fact that the internal absorptions we calculate 



here (see Section 5.7) are, in general, lower than those pre- 
viously used. For example, if we adopt a typical, constant, 
A (Ha) = 1.1 mag value for the internal absorption, the me- 
dian SFR becomes 0.381q;24 Mq yr~^, more consistent with 
other values found in the literature. 

To study the range of SFRs in the different subsamples, 
we binned the data into 0.5 magnitude bins and calculated 
the range (the difference between the maximum and min- 
imum values) for each bin. We find maximum ranges be- 
tween 2.0 and 2.9 dex for the field, SINGS large and Virgo 
subsamples. In every instance, the peak of maximum range 
is located around M(B) = 20.0 magnitudes. On the other 
hand, the range of values in the SINGS small sample, as can 
easily be appreciated from Figure [S] is much lower, with a 
maximum value of 1.3 dex. The same results, for the whole 
sample, are shown in the lower plot of Figure [7] as the nor- 
malized root mean square distance of each data point to the 
best fit. The largest dispersion occurs between M(B) ^ — 17 
and M{B) - -20. 

Figure [9] presents the measured SFR as a function of 
the Hubble stage of each galaxy. We find that most of the 
star forming members of our sample are spiral galaxies with 
stage T ^ 4 (Sb/Sbc). A secondary peak is found for late- 
type spirals (T ^ 7). These results are in good agreement 
with previous studies (see, for example, James et al.||2 004| . 
While there is an expected decline in the SFR toward late- 
type and irregular galaxies the same is not true for very 
early ellipticals. A slight increase in SFR appears to exist 
for early-type galaxies (T ^ —2). This behavior seems to 
contradict the widely accepted theory that elliptical galax- 
ies have exhausted most of their gas reservoirs and do not 
form stars. A series of recent studies (see, e.g., Sarzi et al. 



2006 2010 ) has shown that these galaxies may contain more 



ionized gas than traditionally thought. However, this may be 
mostly due to internal shocks and not necessarily to the star 
formation processes occurring in the galaxy. The number of 
early- type galaxies in our sample is, unfortunately, not large 
enough to make a definitive statement in this debate. 



Figure 10 shows S (the SFR surface density, in units 
of Mq yr~^ kpc~^) of each galaxy plotted against the corre- 
sponding SFR. To compute the size of the galaxy we used 
the optical size D25 from RC3 and converted it to kpc using 
the distance from Table 2. We find that galaxies with higher 
S tend to host larger amounts of star formation, with a 
regression fit of log(5]) = 0.371og(SFR) - 2.28. [KnapeiT^ 



^ The uncertainties are calculated by selecting a random sub- 
sample from the complete set of SFR values and calculating its 
median. We repeated this process a thousand times and selected 



the maximum and minimum median values obtained as the upper 



and lower uncertainties. 
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[Jam es' ("2009) (their Figure 4) published a similar plot for 
the galaxies in the HaGS sample. In their case, however, S 
remains constant for SFRs lower than ^ 1 Mq yr~^. 

As in the case of the SFRs, the result for S shows a 
wide dispersion, with values ranging across three orders of 
magnitude. This can be said for each one of the subgroups 
which make up our sample. 

We also compared our SFRs with the M{K) magnitudes 
of the galaxies in our sample, following the same process as 
explained above. We find that the general behavior of the 
plotted data does not change when using M{K) instead of 
M{B) for the SINGS Large and Smah and for the Virgo 
subsamples. In the case of the Field subsample, the data 
show a higher dispersion when using M(K). 



6.2 Equivalent Avidths 

Figures pT] and p^ show, globally and separated by subsam- 
ple, the EWs of the galaxies in our sample as a function of 
their absolute magnitude in the B-hand. 

We find that our EW values are distributed over almost 
two orders of magnitude, with a median value of 27^20 ^• 
The different subsamples are consistent with this global 
result showing median values that range between 2llg A 
(SINGS Large) and Slt^l^ (SINGS Small). 

The relation between EW and M(B) is basically flat, 
which indicates that, in most cases, an enhanced rate of star 
formation comes with a similar increase in the continuum 
emission produced by evolved stars. This behavior has been 
suspected for almost three decades (see, for example Kenni- 
cutt & Kent 1983) and confirmed in recent studies such as 



Gavazzi et al.| ([1998]) ( t heir F igure 5), [James et al.| ( |2004 ) 
and Knapen fc James] (|2009), the latter two using HaGS 



data. They also detected (see Figure 3 in Knapen & James 
|200 9) a significant number of very luminous galaxies with 
low EW. Their hypothesis that these galaxies are, in their 
majority, massive elliptical and early- type objects with al- 



most no star formation is also supported by Lee et al. ( 2007 ) 



and Lee et al. (2009a), who showed that this result is also 



valid when comparing the EW with the rotation velocity 
of the galaxies. The downturn in EW is not as conspicuous 
in our sample as in these studies, but it seems to exist for 
M{B) < — 17 mag, a region populated mostly by field and 
Virgo galaxies. 

Four out of the ten highest EW values are field galaxies, 
the top one being PGC 43211, a spheroidal galaxy with an 
EW of 880 A. Almost no SINGS galaxies are found among 
the galaxies with the highest EW values. This indicates that, 
while these other galaxies contain SFRs usually above the 
mean for field galaxies, their masses of old stars are pro- 
portionally high. We also find that two of the lowest EW 
measurements, with EWs below 5 A, are intermediate- type 
(Sab and Sb) galaxies from the SINGS large subsample. 
NGC4470, of Sa class and EW - 0.1 A, is the galaxy with 
the lowest EW in our sample. 



Figure 13 (which shows the EW plotted as a function of 
the Hubble stage) does not show any significant correlation 
between EW and morphological type except, maybe, a cer- 
tain increase of the EW values for galaxies of Sb type and 
later. 



7 OTHER STUDIES 

Our current study is one of many surveys based on Ha as a 
tracer of massive stars. In the last decade alone, many sam- 
ples have been observed and reduced using techniques simi- 
lar to those presented here (e.g., [Gavazzi et al.|2002[ Boselli 



eraL][2002| [Gavazzi et al.||2003| [Hunter fc Elmegreen||2004| 



James et a 



et al, 



2005 



. [20041 [Meyer et al 



Gavazzi et al. 



2006 



2004| [Verdes-Montenegro[ 



Kennicutt et al. 



2008) 



Notwithstanding the common goal of characterizing massive 
star formation in the local Universe, these samples have used 
different selection criteria and thus all suffer from different 
biases. Some of the most representative SFR studies in re- 
cent years are addressed below and their results compared 
with some of ours. 

The GOLDMine projecljj has observed and collected 
data for over 3000 galaxies belonging to, among others, the 
Coma supercluser and the Virgo, Cancer and Hercules clus- 
ters. This data includes Hi, IR and optical continuum mea- 
surements, as well as spectral observations. Two large sets 
of original Ha data have been published by Gavazzi and col- 
laborators. Gavazzi et al. (2002|) observed 369 galaxies with 



types later than SOa in the Virgo and Coma/A1367 clus- 
ters. iG^^zzT^r^ (2006) extended the sample with some 
273 galaxies in the Virgo and Cancer clusters and in the 
Coma/A1367 and Hercules superclusters. 



We have already introduced the SINGS survey ( Ken- 
nicutt et al.[[2003) in Section [2] Aimed to characterize the 



infrared emission in the local Universe, it includes 75 galax- 
ies within 30 Mpc from elliptical to irregular and with a 
large range of SFRs. The SINGS team has obtained Ha and 
continuum broad-band imaging for the whole sample but, 
although the images have been made public, no catalogue of 
galaxy-integrated SFRs or EWs has been published. 

Our survey has also made extensive use of the HaGS 
work (James et al."2004), which observed 334 galaxies ex- 
tracted from the Uppsala Galaxy Catalogue (Nilson 1973} 
of types later than SOa and with m(B) < 14.5 in the Palo- 
mar observatory plates. The targets were selected to include 
galaxies with recession velocities from to 3000 km s~^. As 
in our case, they found a clear trend between SFRs and ab- 
solute magnitudes and an apparent lack of correlation of EW 



with luminosity and morphological type (see also Knapen & 
James[[2009| ). 



Although it is generally accepted that most of the mas- 
sive star formation in galaxies occurs in the spiral arms 
of late- type galaxies, some authors have questioned the 
necessity of a spiral density wave to achieve an impor- 
tant amount of star formation (see, for example, [Hunter, 



Gallagher Sz Rautenkra nz[|1982[ [ Elmegreen fc Elmegreen 
1986 



Elmegreen 2011 and references therein). Hunter & 
Elmegreen ( ,2004 hereafter HE04) focused on observing a 
sample of 142 mainly irregular galaxies. Some Sm and blue 
compact dwarfs (BCD) were observed for comparison. The 
sample of irregular galaxies was compiled from the |Fisher| 
& Tully ( 1975| ) Hi atlas of galaxies and thus includes only 
galaxies containing gas. They found Ha emission in over 80 
of their 96 irregular galaxies and in all the Sm and dwarf 
(BCD) galaxies except one. The SFRs are low, with a mean 
value of 0.1 Mq yr~^ and a maximum of 0.5 Mq yr~^. 



^ http://goldmine.mib.infn.it/ 
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The Survey for Ionization in Neutral-Gas Galaxies 
(SINGG) aims to observe 468 galaxies selected from the 
Hi Parkes All Sky Survey (HIPASS, p'eyer et al.|2004t with 
peak flux densities above 0.05 Jy. The galaxies have been 
selected to cover a range of log(MHi/M0) from 8 to 10.5. 
The most recent data release we have been able to find is the 
one published by HanishetaL ( ^200 6 ) which contains obser- 
vations and derived values for 111 galaxies. Among those, 
every galaxy with Mhi > 3 x 10^ Mq has also been detected 
in Ha, thus constraining the amount of gas necessary to 
trigger global star formation in a galaxy. 



Kennicutt et al.| f2QQ8) (hereafter RCK08) have made, 
along with HaGS, one of the first attempts to observe a 
volume-limited sample of galaxies using Ha. Two main fac- 
tors make this a difficult task: firstly, the mapping of galaxies 
(especially irregular ones) in the local universe is currently 
not complete; and secondly, even for the known galaxies, the 
high uncertainties in the distances make it difficult to define 
the set of galaxies within a certain radius. While these prob- 
lems cannot be solved at the present moment, RCK08 man- 
aged to observe a sample of 261 S and Irr galaxies within 
11 Mpc with galactic latitudes above 20° and ui{B) ^ 15 
magnitudes, and an additional set of 175 galaxies with lower 
galactic latitudes or fainter. This sample is shared with other 
legacy surveys to be carried out using Galex and Spitzer 
which will provide interesting complementary information 
in the near and far UV and in the IR (3.6 to 160 /xm). 

Finally, 'Karachentsev fc Kaisin| ( |2010 ) presented a 
study of 52 galaxies observed in Ha as part of a larger sur- 
vey aimed to study a sample of galaxies in the Local Volume 
(207 galaxies within 10 Mpc). An important number of the 
galaxies observed belong to the Maffei 2/IC 342 complex. 
The sample contains galaxies of all morphological types but 
with prevalence of irregular galaxies. 

Figure [14] shows a series of histograms comparing the 
number of galaxies versus the morphological class (in the 
form of Hubble stages) for our sample and the different 
samples presented above. Most of these studies lack a good 
representation (or any representation at all) of early-type 
galaxies. While this is a natural thing in surveys like HE04, 
aimed to study irregular galaxies specifically, it indicates 
that a general bias towards late-type (SO and later) galax- 
ies dominates the literature of Ha surveys. This comparison 
also shows that H i-selected samples, as SINGG and our own 
sample, span a wider range of morphological classes. 

Figures [15] to [T7| present similar histograms for distance, 
SFR and EW respectively. SINGS does not provide measure- 
ments of Ha fiuxes (and, consequently, L(Ha) and SFR) or 
EW. EWs for HE04 and RCK08 are not available in the 
public catalogues. The EWs from RCK08 we use to com- 



plete our data (see Section 5.7) have been taken directly 



from the headers of the images, but they are not included 
in the catalogue published in CDS. 

The distance histograms show that most of the galax- 
ies observed to date with resolved Ha imaging are closer 
than 50 Mpc even when, among the surveys discussed here, 
only RCK08 and HaGS include the distance as a selection 
criterion. The exceptions are SINGG, which observed galax- 
ies as far as 



the clusters observed (^ lOOMpajfor the Coma Cluster, 
16.4 to 22.9 Mpc for the Virgo Cluster - [Mei et al.|20Q7l and 
> 150 Mpc^ for the Hercules Cluster). The galaxies in our 
sample have a mean distance of 15.8 Mpc, similar to those 
in the rest of the surveys, except for RCK08 with 6.8 Mpc 
and [Gavazzi et al." (2002", 2006) with - 50 Mpc. 

All the SFR histograms shown in Figure [16] tell a sim- 
ilar story, with most galaxies hosting small amounts of star 
formation (less than O.5M0yr~^). The number of galax- 
ies with rates of star formation higher than ^ lM0yr~^ 
drops dramatically in every case. [Gavazzi et aT] ( |2002| show 
a significant number of galaxies with high SFRs (up to 
^ 6 Mq yr""*^). This may be due to a systematic overestima- 
tion of their Ha fiuxes, since the distribution of morphologi- 



cal types (see Figure 14) does not indicate an overwhelming 
presence of highly efficient star forming galaxies. Eleven of 
the SINGG galaxies present SFRs higher than lOMoyr"^ 
with the extreme case of M 83 with a SFR of 145 M© yr~\ 
Those cases are suspect, as such behaviour is not found in 



any of the other surveys. Hunter & Elmegreen (2004) and 



Karachentsev & Kaisin ( 2010 ) show virtually no SFRs above 
0.5 M0 yr~ , which can be attributed to the fact that their 
samples are mostly or totally composed of irregular galaxies 
with very low SFRs. 



A similar behaviour is found (Figure 17) for the EW of 
the SINGG galaxies. The large number of galaxies with very 
high EW (> 140 A, with a maximum value of 450 A for the 
peculiar galaxy NGC 5253) is surprising. These exceptions 
notwithstanding, the comparison between the different his- 
tograms show a similar pattern, with mean EWs between 20 
and 35 A (in the case of SINGG the mean is 50.7 A but this 
is mainly due to some extremely high values). 

Figure [18] shows the comparison between our SFR data 
and the values found in the literature. Except for a global off- 
set, which can be explained by the different internal absorp- 



tion used in this work (see Section 6.1 ), the global agreement 
is good, with a fit of slope 0.90 and a correlation coefficient 
R^ = 0.82. This good agreement is also valid if we compare 
our data with each individual study. We find a slope of 0.95 
and R^ = 0.98 for HaGS, 1.01 and 0.88 for RCK08, 1.35 and 
0.92 for CM and 0.86 and 0.58 for SINGS. The excellent fit 
in the case of HaGS is especially important as our aperture 
selection and fiux measurement methods are essentially the 
same. As we have already commented, the SINGS team has 
not published SFR measurements based on their Ha images. 



Kennicutt et al.|(|2003) gives the SFR values from the IRAS 
catalogue (Fullmer & Lonsdale 1989), based not on Ha ob- 



servations but on IR measurements from 60 to 100 /im. This 
explains that the fit using SINGS values is worse than those 
using the other points. In any case, these comparisons show 
that the differences between our SFR values and those in 
the literature lie well within the typical errors for this kind 
of measurement, whatever the methods applied and the spe- 
cific luminosity-to-SFR constant used in each study. 



Mpc, and Gavazzi et al. 



most 
of whose galaxies are located at the specific distances of 



(|2002^^2006) 



^ From NED 
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8 CONCLUSIONS AND FUTURE WORK 

We have presented new Ha and i?-band observations of 72 
out of the 156 galaxies sampled in the JCMT Nearby Galax- 
ies Legacy Survey. This sample is partly Hi-selected and 
aims to include galaxies within 25Mpc with morphological 
types from E to Irr. The sample is largely formed by field 
galaxies but also includes targets from the SINGS survey 
and a fair amount of Virgo cluster galaxies. 

For the remaining 84 galaxies we have compiled images 
from the literature. We processed all these data using well- 
defined methods and criteria and present Ha fluxes, SFR 
and EW measurements for the whole sample. With this ar- 
ticle, we make our new data publicly available. 

We analyzed the data reaching, among others, the fol- 
lowing conclusions: 

(i) The SFR shows a clear increasing trend when plotted 
against M(B). This trend is dependent on the subsample 
considered, with a better fit for the field and SINGS Small 
galaxies than for those in the SINGS Large or Virgo sub- 
samples. 

(ii) We find a median SFR across our sample of 
0.2 Mq yr~^, while in the case of the galaxies in the SINGS 
Large subsample the median SFR is 1 Mq yr~^. 
(iii) The SFR values range across almost three orders of 
magnitude, except in the case of the SINGS Small galax- 
ies where the values are confined to a range of one order of 
magnitude. 

(iv) The SFR of our sample peaks for Sb/Sb-c galaxies (T ^ 
4) with a secondary peak for late- type Sc galaxies (T ^ 7). 

(v) We plotted EW versus M{B) confirming that the rela- 
tion between EW and luminosity is practically fiat. This is 
also valid when comparing the EWs with the morphologi- 
cal type, although a trend seems to exist leading to slightly 
higher EW for late- type and irregular galaxies. 
(vi) The comparison of our SFR values with those found in 
the literature shows an excellent agreement for SFR, well 
within the typical margins of error. 

This paper has presented our Ha line observations, as 
well as a general study of the star forming properties of 
the sample and a comparison with previous similar studies. 
However, this is only the first step towards the more am- 
bitious goal of comparing the star formation in the galax- 
ies of the sample with their gas content, both atomic and 
molecular. ^^C0J = 3 — 2 maps of the whole sample are 
available as a part of the NGLS observations carried out 
with HARP-B in the JCMT, thus allowing a global com- 
parison between L(Ha) and L(CO) for all our galaxies. This 
will be the largest study correlating the star formation with 
the amount of warm, dense molecular hydrogen present in 
a galaxy and may be key to understanding the physics in- 
volved in large-scale star formation in the nearby universe. 

Furthermore, the rich variety of literature data avail- 
able for a large part of our sample, from IR data obtained 
with Spitzer to FUV from GALEX as well as Hi from a 
variety of sources ( THINGS: [Walter et al 



van der HuIst]|2QQ2 



2008 



^^ ^ WHISP: 

etc.) and CO J = 1-0 and J = 2 - 1 
(e.g., HERACLES: (Leroy et al.|2009t allows us to study the 
behavior of the different gas and dust components over the 
complete morphological range. 

Finally, the JCMT NGLS has been awarded guaranteed 



time to observe the whole sample using SCUBA-2 to obtain 
continuum imaging in 450 and 850 /xm and study the cold 
dust within galaxies with an unprecedented sensitivity and 
resolution. These data will also be compared to the current 
Ha and other data in forthcoming papers. 
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Figure 1. Histogram of the number of galaxies as a function of the Hubble stage, T, as listed in the Third Reference Catalogue of Bright 
Galaxies (RC3). The four different subsamples are shown with different patterns and colors. 
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Figure 2. Ha continuum-subtracted images of the field galaxies of our sample. The ellipse shows the aperture used to perform the 
photometry (see Section [5^. All the images are oriented North up, East to the left. The right bottom line is always SOarcsec. 
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Figure 2. Continued 
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Figure 2. Continued 
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Figure 2. Continued 
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Figure 2. Continued 
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Figure 2. Continued 
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Figure 3. As Figure |2] but now for the SINGS Large galaxies in our sample. 
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Figure 3. Continued 
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Figure 3. Continued 
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Figure 4. As Figure [2] but now for the SINGS Small galaxies in our sample. 
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Figure 4. Continued 
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Figure 5. As Figure [2] but now for the Virgo Small galaxies in our sample. 
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Figure 5. Continued 
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Figure 6. Relation between the transmission of the Ho; and Nil hnes for different narrow filters used in this work, as a function of the 
filter width. The dashed line indicates the exponential fit obtained (see Section [5^. 
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Figure 7. Logarithm of the SFR (in units of solar masses per year) versus the absolute B magnitude of the galaxies in our sample. The 
four subsamples which constitute the complete list are marked with different colors and symbols. The black solid line indicates the least 
squares regression fit for the plotted data. The lower plot shows the normalized root mean square of each galaxy to its fitted value. 
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Figure 8. As FigurelTlbut now with the subsamples divided into different panels. The dashed lines indicate the mean values of log(SFR) 
for each subsample while the dotted ones indicate the median values. 
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Figure 9. Logarithm of the SFR (in units of solar masses per year) versus the Hubble stage of the galaxies in the sample. The different 
subsamples in our complete list are plotted separately. The solid black line shows the mean values of the data binned to AT = 1. The 
dashed blue line shows the median values calculated with the same procedure. 
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Figure 10. As Figure [7| but now plotting the logarithm of E in units of solar masses per year and per square kiloparsec versus the SFR 
in units of solar masses per year. 
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Figure 11. As Figure [t] but now with the logarithm of the equivalent width in Angstrom. The dashed line indicates the mean value of 
log(EW) for the full sample while the dotted one indicates the median value. 
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Figure 12. As Figure [TT] but now with the subsamples divided into different panels. 
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Figure 13. As Figure |9] but now with the logarithm of the equivalent width in Angstrom. 
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Figure 14. Histograms comparing the RC3 Hubble stage distribution, T, of our sample (top left panel) with several samples found in 
the literature: the [Gavazzi et aL] ( |2002||2006|> survey of nea rby galaxy clusters, the S INGS sample ([Kennicutt et al.|2003| , the sample of 
irregular and Sm galaxies observed by [Hunter fc Elmegreen| ( |2004|), the HaGS survey ( [James et al.[2004[ ), the SINGG survey ([Meurer et aP 
2006|>, the volume-limited sample within 11 Mpc studied by [Kennicutt et aL] ( [2Q08] > and the 52 galaxies observed in Ha by [Karachentsev 
Sz Kaisiii] ([2010 ) as part of a larger Local Volume survey. Note that the last bin in each histogram includes the galaxies with values of T 
larger than the upper limit in the plot. Further details are given in the text of Section [7| 
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Figure 15. As Figure [14] but now showing the histograms of the distance in units of Mpc. 
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Figure 16. As FigurelTSbut now showing the histograms of the SFR in units of Mq yr ^ . The values plotted here have been calculated 
directly from the F{Yia) measurements taken from the catalogues using the same extinction and F(HQ;)-to-SFR values as in Sections |5.5| 
andO 



© 2011 RAS, MNRAS 000, [Tp5 



46 J. R. Sanchez- Gallego et al. 



30 
25 
20 
15 
10 
5 


70 
60 
50 
40 
30 
20 
10 


80 
70 
60 
50 
40 
30 
20 
10 




i. 


f 


I 










NGLS 


- 


i 


1Z 

1. 


— 




\ 




3 


3 




r^ 



n \ 1 r 



J 



n r- 

HaGS 



1 1 \ r 



RCK08 



^^ 



n \ r 




20 40 60 80 100 120 140 

EW [A] 



20 40 60 80 100 120 140 

EW [A] 



Figure 17. As Figure |l4| but now showing the histograms of the EW in Angstroms. 
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Figure 18. Comparison of SFR values calculated in this work with those obtained from the literature. Except in the case of the SINGS 
values, all the SFRs have been derived from the same images used in the literature. The solid line indicates y = x. Further details can 
be found in the text. 
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